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CZ3 ■ Abstract 

O 

c/5 ■ In this work we suggest a simple modification of tlie Newton's cooling law that can model 

^~»! Mpemba effect. We introduce, in the usual Newton's law, i.e. linear differential equation, 

o ' an additional term proportional to the quadrate of the geometrical average value of the ini- 

tial and latter difference between liquid and cooling thermostat (environment) temperature. 
It, after simple transformations, yields usual Newton's linear differential equation but with 
!>■ I modified cooling parameter. This modified cooling parameter represents sum of the usual 

cooling parameter and an additional term directly proportional to the difference between 
initial temperature of the liquid and cooling thermostat temperature. Corresponding solu- 
tion of the modified Newton's cooling equation, i.e. temperature decrease during time, has 
lO ' an additional exponential term with negative argument proportional to mentioned difference 

between initial temperature of the liquid and temperature of the cooling thermostat. (It 
can be observed that appearance of such "hysteresis" or "memory" terms can imply some 
circulation effects by liquid flow but it goes over basic intention of this work.) It admits that 
liquid with higher initial temperature, after a characteristic time interval, cools faster than 
the same liquid with lower initial temperature, in full agreement with experimental data on 



C^ ' the Mpemba effect. 



In this work we shall suggest a simple modification of the Newton's cooling law [1] that can 
model Mpemba effect [2]- [5]. We shall introduce, in the usual Newton's law, i.e. linear differential 
equation, an additional term proportional to the quadrate of the geometrical average value of the 
initial and latter difference between liquid and cooling thermostat (environment) temperature. It, 
after simple transformations, yields usual Newton's linear differential equation but with modified 
cooling parameter. This modified cooling parameter represents sum of the usual cooling parame- 
ter and an additional term directly proportional to the difference between initial temperature of 
the liquid and cooling thermostat temperature. Corresponding solution of the modified Newton's 
cooling equation, i.e. temperature decrease during time, has an additional exponential term with 



negative argument proportional to mentioned difference between initial temperature of the liquid 
and temperature of the cooling thermostat. (It can be observed that appearance of such "hystere- 
sis" or "memory" terms can imply some circulation effects by liquid flow but it goes over basic 
intention of this work.) It admits that liquid with higher initial temperature, after a characteristic 
time interval, cools faster than the same liquid with lower initial temperature, in full agreement 
with experimental data on the Mpemba effect. 

As it is well-known [1] usual Newton's law of cooling represents the following simple linear 
differential equation 

^=-a(r-T,) (1) 

with simple solution 

T = re + (To-Te)exp[-at]. (2) 

Here T represents the temperature of a liquid with initial value Tq, Tg - constant, smaller than 
T, temperature of a cooling wall (thermostat or environment) in the contact with liquid, while a 
represents a cooling parameter that simplifiedly expresses thermodynamics of the cooling process. 
It is supposed implicitly that cooling parameter is independent of the initial temperature of the 
liquid as well as of the cooling wall temperature. 

Consider now described by (1) or (2) thermodynamical interaction between liquid and cooling 
wall in two different cases, precisely for two different initial temperatures of the liquid, first one 
Tio and second one T20 so that Tio is larger than T20 , i.e. 

Tio > T20. (3) 

Then (2), (3) imply 

Ti — Te _ Tio — Te , . 

T2-T,~ T20 - Te ^ ^ 

where Ti and T2 represents the temperature of the liquid in the first case and second case in some 

(finite) later time moment t. It means that temperature of the liquid in the first case will be larger 

than the temperature of the liquid in the second case in any (finite) later time moment. 

Previous conclusion implies unambiguously that usual Newton's cooling law (1) cannot model 
experimental data of the Mpemba effect [2]- [5] according to which, simply speaking, hot liquid 
freezes faster than cold liquid. 

Also, it is well known [1] that there are many different non-linear modifications of the Newton 
cooling law. For example there is the following modification 

f = -<.(T-r.r (5) 

for 1.3 < n < 1.6 in the Taylor' s model and for n = | in the Dulong- Petit model. Or, there is 
so-called Newton - Stefan modification of the usual Newton's cooling law 

^ = _a(T - Te) - h{T - T,f - c{T - T,f - d{T - T,)\ (6) 

where cooling parameters a, h and c are dependent of Tg. 

We shall suggest a simple, quasi-quadratic modification of the usual Newton's cooling law 

^ = _a(r-Te)-6<T-Te>2 (7) 



where 

<T-T,>= ((T - Te)(To - Te))^ (8) 

represents the geometric average value of the initial and latter difference between liquid and cooling 
wall temperature. 

Introduction of (8) in (7) yields 

^ = -(a+6(To-T,))(T-T,) (9) 

representing a linear differential equation similar to usual Newton's cooling law (1). However, in 
distinction to usual Newton's cooling law, (9) holds cooling parameter (a + 6(To — Tg)) dependent 
of the initial difference between liquid and cooling wall. 

As it is not hard to see equation (9) holds the simple solution 

T = Te + (To - Te) exp[-(a + b{To - T,))t]. (10) 

It can be observed that in the limit of small (Tq — Tg), i.e. for (To^^e) ^ f , modified Newton's 
cooling law, i.e. equation (9) and its solution (10) turn out usual Newton's cooling law, i.e. in the 
equation (1) and its solution (2). 

Consider now described by (9) or (10) thermodynamical interaction between liquid and cooling 
wall in two different cases, precisely for two different initial temperatures of the liquid, first one 
Tio and second one T20 so that Tiq is larger than T20 , i.e. 

Tio > T20. (11) 

Then, according to (9)-(ll), between temperature of the liquid in the first case Ti and tem- 
perature of the liquid in the second case T2, there is the following relation 

^^ = ^^exp[-6(Tio -T2o)t]. (12) 

It represents a time decreasing function that is, according to (11), initially larger than 1 but that 
after a characteristic time moment 

T = —^ tM^^^^^] (13) 

becomes smaller than 1. In other words, after characteristic time moment r, Ti becomes smaller 
than T2 even if TiqT is larger than T20. 

As it is not hard to see, all this can be used as a simple model for the Mpemba effect [2]- [5]. 

It can be observed that average temperature difference < T — Te > (8) representing geometrical 
average value of the initial and later temperature difference does not represent unique form of the 
average temperature difference that introduced in (7) can model the Mpemba effect. For example, 
as it is not hard to see (but which will not be analyzed here explicitly), arithmetical average value 
of the initial and later temperature difference introduced in (7) can model the Mpemba effect 
too. But, in this case, expression (7) represents the differential equation with constant, linear and 
really existing quadratic term depending of T — Tg as the variable and Tq — T^. as the parameter 
(which implies that solution of such equation decreases during time faster than (10) ). 

In conclusion it can be repeated and pointed out the following. In this work we suggest a 
simple modification of the Newton's cooling law that can model Mpemba effect. We introduce. 



in the usual Newton's law, i.e. linear differential equation, an additional term proportional to 
the quadrate of the geometrical average value of the initial and latter difference between liquid 
and cooling thermostat (environment) temperature. It, after simple transformations, yields usual 
Newton's linear differential equation but with modified cooling parameter. This modified cool- 
ing parameter represents sum of the usual cooling parameter and an additional term directly 
proportional to the difference between initial temperature of the liquid and cooling thermostat 
temperature. Corresponding solution of the modified Newton's cooling equation, i.e. temperature 
decrease during time, has an additional exponential term with negative argument proportional 
to mentioned difference between initial temperature of the liquid and temperature of the cooling 
thermostat. (It can be observed that appearance of such "hysteresis" or "memory" terms can 
imply some circulation effects by liquid flow but it goes over basic intention of this work.) It 
admits that liquid with higher initial temperature, after a characteristic time interval, cools faster 
than the same liquid with lower initial temperature, in full agreement with experimental data on 
the Mpemba effect. 
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